Dynamics of spin-1 bosons in an optical lattice: spin mixing, quantum phase revival 
spectroscopy and effective three-body interactions 
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We study the dynamics of spin-1 atoms in a periodic optical-lattice potential and an external 
magnetic field in a quantum quench scenario where we start from a superfluid ground state in 
a shallow lattice potential and suddenly raise the lattice depth. The time evolution of the non- 
equilibrium state, thus created, shows collective collapse-and-revival oscillations of matter-wave 
coherence as well as oscillations in the spin populations. We show that the complex pattern of 
these two types of oscillations reveals details about the superfluid and magnetic properties of the 
initial many-body ground state. Furthermore, we show that the strengths of the spin-dependent and 
spin-independent atom-atom interactions can be deduced from the observations. The Hamiltonian 
that describes the physics of the final deep lattice not only contains two-body interactions but also 
effective multi-body interactions, which arise due to virtual excitations to higher bands. We derive 
these effective spin-dependent three-body interaction parameters for spin-1 atoms and describe how 
spin-mixing is affected. Spinor atoms are unique in the sense that multi-body interactions are 
directly evident in the in-situ number densities in addition to the momentum distributions. We 
treat both antiferromagnetic (e.g. 23 Na atoms) and ferromagnetic (e.g. 87 Rb and 41 K) condensates. 

PACS numbers: 03.75. Mn, 03.75.Kk, 67.85.-d, 67.85.Fg 



I. INTRODUCTION 

Quantum degenerate ultracold atoms with spin-degree 
of freedom exhibit both magnetic order and superfluidity, 
offering a rich system in which to explore quantum coher- 
ence, long-range order, magnetism and symmetry break- 
ing. Many aspects of spinor atoms in a trap have been 
investigated with spin F = 1 atoms, such as 23 Na and 
87 Rb (Ml. Spin-2 @-|l and spin-3 @, []jj spinor gases 
have been studied to a lesser extent. Spinor condensates 
are described by a vector order parameter [Ill4l3| . The 
distinctive feature is its spin-dependent interaction which 
organizes spins giving rise to ferromagnetic and antifer- 
romagnetic (polar) order. It can also coherently convert 
a spin m = 1 and a m = — 1 atom to two m = atoms 
and vice versa [l44l6j . while conserving magnetization 
and energy. 

In parallel to ultracold spinor physics, optical lat- 
tices have become a powerful tool to create strongly cor- 
related many-body states of bosons and fermions fl7r - 
l2l| . Lattice systems offer flexibility as the lattice pa- 
rameters and particle interactions can be controlled eas- 
ily, thereby facilitati ng pro gress towards the creation of 
quantum emulators [221 , l23l ] . Since the seminal observa- 
tion of the superfluid to Mott insulator transition with 
spinless bosons [I5], steady progress is being made to- 
wards the understanding of spinor atoms in an optical 
lattice [24426J ■ Issues of temperature and entropy [27| are 
among the challenges that need to be overcome to create 
a many-body correlated state of spin-1 atoms. Theo- 
retical studies of lattice-trapped spinor condensates have 
mainly explored the phase diagram and the nature of the 
superfluid- Mott insulator transition |28l - [33| . 

Due to the tunability of cold atom and optical lattice 
parameters, it is also possible to study non-equilibrium 



dynamics. Dynamics of many-body quantum systems is 
still an emerging field, and only a number of issues have 
so far been investigated [34| . An early experiment [35| 
studied the dynamics of spinless bosons in a suddenly- 
raised optical lattice, observing the collapse and revival 
of the matter wave field in the momentum distribution. 
In a more-recent experiment [361 ] , tens of oscillations in 
the momentum distribution or visibility were observed, 
and the predicted [33] signature of effective higher-body 
interactions confirmed. As for spin-1 atoms, dynamical 
studies have mainly focused on large atom continuum or 
trapped systems in the mean field regime exploring spin- 
mixing dynamics [14 , Il5| , quantum quench dynamics [38l . 
l39j , and various instabilities J4Q] . 

The goal of this paper is to study the dynamics of 
spin-1 bosonic atoms in a three-dimensional (3D) optical 
lattice and probe its many-body state and system prop- 
erties. Starting with a ferromagnetic ( 87 Rb) or antifer- 
romagnetic ( 23 Na) superfluid ground state in a shallow 
lattice, suddenly raising the lattice depth creates a non- 
equilibrium state, which can be followed in various sce- 
narios - with and without a magnetic field and with and 
without effective three-body interactions. The evolution 
shows collapse and revival of matter- wave coherence mea- 
sured by visibility oscillations, in a more complex pattern 
than for spinless bosons [35|, EH • It also shows oscillations 
in spin populations due to the combined effect of the spin- 
mixing collisions of the m = and m = ±1 components 
and differential level shifts proportional to the square of 
the magnetic field strength, the quadratic Zeeman shift. 
Linear Zeeman shifts do not affect the behavior of spinor 
condensates. Both spin-mixing and visibility oscillations 
reveal details about the system such as the composition 
of the initial many-body state, and thereby its superfluid 
and magnetic properties. 



By analyzing the frequency spectrum of the visibil- 
ity, we show that the ratio U 2 /Uq of spin-dependent 
and spin-independent atom-atom interactions can be de- 
duced. Combined with spectra of spin-mixing dynam- 
ics at various magnetic field strengths, this gives us a 
method to measure the interaction couplings for spin-1 
atoms. Finally, we find that the presence of quadratic 
Zeeman shift enhances spin mixing oscillations for ferro- 
magnets and shows collapse and partial revivals in the 
transverse magnetization. 

The Hamiltonian that quantitatively describes the 
physics of the final deep lattice comprises of two-body 
as well as effective multi-body interactions, which arise 
due to virtual excitations to higher bands. We derive 
the induced three-body interaction parameters for spin-1 
atoms in a deep harmonic well, approximating the mini- 
mum of a single lattice site as such, and find the existence 
of spin- dependent three-body interactions. We show how 
to detect the signature of the three-body interactions and 
argue that they arc directly exhibited in the in situ den- 
sity as opposed to the time of flight visibility measure- 
ments as is the case for spinless bosons [36| . 

The article is organized as follows. In Sec. [II] we setup 
the spin-1 Bosc-Hubbard model, sketch the mean-field 
theory to determine the initial ground state, describe the 
exact Hamiltonian after the quench, and discuss observ- 
ables and computational aspects. We present our main 
results in Sec lIIIlUVl andlVl Section lHll explores the non- 
equilibrium dynamics of antifcrromagnctic spin-1 atoms, 
with and without a magnetic field. Section [IVI describes 
the dynamics for a ferromagnetic spinor. ScctionlVlshows 
the effects of effective three-body interactions in the dy- 
namics. We summarize our results in Sec. IVI1 A deriva- 
tion of the effective three-body interaction is given in the 
appendix. 



II. MODEL AND COMPUTATIONAL ASPECTS 
A. Shallow Lattice Hamiltonian 

Ultracold spin-1 bosons in the lowest band of a 3D 
cubic optical lattice and an external magnetic field B 
along the z axis are modeled by the free energy 
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Here a\ m is the creation operator of a boson in magnetic 
sublevel m = — 1, 0, or 1 in the energetically- lowest Wan- 
nier function or orbital of lattice site i. The first term in 
Eq. ([TJ represents the hopping of atoms between nearest 
neighbor sites and is proportional to the hopping energy 



J. The second term describes the on-site atom-atom re- 
pulsion with strength Uq > 0, hi = J2 m him, and n im = 
a i m a im- The third term is a spin-dependent atom- atom 
interaction with a strength U 2 that can be either positive 
or negative. The three operators Ti = (Txiifyiifzi) on 
site i satisfy angular momentum commutation rules and 
are defined by F ql = J2 m , m ' a l m (F q ) mm >a im , for q = x, 
y, or z and (F q ) mm i are matrix elements m, ml of compo- 
nent a of the spin-1 angular momentum F. The fourth 
term corresponds to the quadratic Zeeman energy of the 
magnetic sublcvcls with strength 5. Finally, the terms 
containing the Lagrange multipliers [i and Hm control the 
total atom number, N = J^ hi, and total magnetization, 
M = J^thii, respectively. Here m, = F zl = J2 m m him 
is the on-site magnetization. Both N and M commute 
with H. (The large linear Zeeman Hamiltonian of the 
atoms is "absorbed" in the term — (J-mM and seen not to 
affect the physics of the spinor condensate.) 

The interaction strengths arc given by Uq = 
47rh 2 n(a + 2a 2 )/(3M a ) and U 2 = 4nft 2 n(a 2 - a )/(3M o ) 
[llj , where a^ with F = or 2 are scattering lengths for 
the two allowed values of the total angular momentum 
of s-wave collisions of two spin-1 particles at zero colli- 
sion energy and zero magnetic field. S'-wave scattering 
with total angular momentum F — 1 is prohibited due to 
bosonic wave function symmetry. The mean density of 
the local orbital n is determined by the laser parameters 
and polarizability of the atom. Finally, M a is the mass 
of the atom, ft, = h/(2ir), and h is Planck's constant. 

The ratio of the two interaction strengths is indepen- 
dent of lattice parameters as the n dependence cancels. 
For 23 Na U 2 /U = +0.036(3) [H and the system is an- 
tiferromagnetic. For 87 Rb U 2 /U = -0.0046(7) [H| and 
the system is ferromagnetic. The quoted uncertainty in 
U 2 /Uq for 23 Na and 87 Rb are one-standard deviation and 
obtained from the corresponding references. For 41 K 
we find U 2 /U = -0.012 0, |H|. For our investiga- 
tion we use U 2 /Uo = 0.04 for sodium and —0.005 for 
rubidium. The quadratic Zeeman strength S = 5qB 2 
with S /h = 27.68 Hz/(mT) 2 for 23 Na and S /h = 7.189 
Hz/(mT) 2 for 87 Rb. 

The phase diagram for the spin-1 Bose-Hubbard model 
has been calculated with numerical methods such as 
Quantum Monte Carlo [3l| and density-matrix renormal- 
ization group [46| for a one-dimensional lattice. Mean- 
field approaches for spin-1 bosons, which give pre dictions 
for any dimension, have also been performed 132.1471 148| 
and are extensions of those for scalar bosons [49|, [5(| ■ As 
mean-field models are most predictive in three dimen- 
sions and we focus on such a lattice, we use the decou- 
pling mean-field theory [32] to find the initial many-body 
ground state. 

In a mean-field approximation, the hopping term in 
Eq. {1} can be decoupled as 
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when fluctuations around the equilibrium value are neg- 



ligiblc. Wc can define ip m = (a,j m ), for m = 1,0,-1, as 
the site-independent supcrfiuid order parameter. Using 
Eq. ([2]) we can rewrite Eq. ([T]) as a sum of independent 
single site Hamiltonians, H = J^, -ff' nf where 
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and z is the number of nearest neighbors, e.g. z = 6 in 
3D. For a given [i and \im the superfluid order parame- 
ters and ground state wavefunction are obtained by find- 
ing those values of tp m for which the energetically-lowest 
eigenstate of Hf 1 is smallest. 

The character of the ground state depends on whether 



the spin-dependent term XJi is positive or negative [11 
For the antiferromagnetic C/ 2 > superfluid ground 
states, the order parameters can be written as ip m = 
y/p^e l9 Dl n0 (a,(3,"/), while for the ferromagnetic U 2 < 
supcrfiuid we have tp m = ^/p^e l6 Dl nl (a,f3,j). Here the 
functions D J MM , (a, /3, 7) are Wigner rotation matrices 
[5l| with Eulcr angles a, (3, and 7 determined by min- 
imizing Hf 1 , The real valued p s and angle 9 arc the 
spin-independent supcrfiuid density and a global phase, 
respectively. We have p s < (h) . 

Within mean-field theory the many-body superfluid 
wavefunction is given by the product wavefunction 
Y[ t \GS)i over sites i, where |GS) = J2a c n\™) anc -^ ^cts 
\n) = |n_i,no,rii) are elements of the occupation-number 
basis of Fock states of the three m projections. The 
single-site wavefunction is a superposition of Fock states 
with amplitudes c r -j. In fact, it is also a superposition of 
Fock states with different magnetization. 

We only present results for ground states with zero 
magnetization (m*) = at every site, ensured by set- 
ting [Am = 0. Results for other magnetizations show 
similar physics. For U 2 > the bosons condense into 
a state with (Fi) = 0. This is called a polar (anti- 
ferromagnetic) supcrfiuid. There are two kinds of po- 
lar order [32j: the ground state is a transverse polar 
state with (tp^i,tpo,ip-i) = y/Jh(l, 0, 1)/V2 when 5 < 
0, and a longitudinal polar state with ("0-1 7 V'OiV'i ) = 
v /p7(0,l,0) when 5 > 0. For ferromagnetic atoms 
with U2 < the magnetic order maximizes the to- 
tal angular momentum with (J r i) 2 = 1 [ll| . leading 
to a partially-magnetized superfluid with order param- 
eters (V>_i,Vo,V'i) = v/p7(l,\/2,l)/2 for < S < 2U 2 , 
and a longitudinal superfluid with order parameters 
^7(0,1,0) for 8 >2U 2 . 

Our numerical simulations are performed in the occu- 
pation number basis. Only basis functions with n_i + 
n + n\ < n max are included. We use n max = 6 leading 
to 84 basis functions in a site and negligible truncation 
errors when the mean atom number per site is less than 



three. All current optical-lattice experiments use mean 
atom numbers of this order of magnitude. 



B. Deep Lattice Hamiltonian 

After preparing the initial superfluid the depth of the 
optical lattice is suddenly increased so that tunneling is 
turned off. This lattice ramp-up is assumed to be slow 
enough to prevent excitations to a higher band yet fast 
enough compared to interactions. We can then treat sub- 
sequent time evolution due to the single-site Hamiltonian 
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exactly. As each site evolves under the same Hamil- 
tonian, we have suppressed the site index. For in- 
duced three-body interactions, additional terms appear 
in H Rnal as discussed in Sec. |V] Recent observation of 
multi-body effects for lattice-trapped spinless bosons [36| , 
where a similar quench was used, a mean-field treat- 
ment [23, [52( of the initial state followed by exact on-site 
evolution was found to agree well with the experiment. If 
lattice sites are not completely decoupled (J^0) during 
the evolution, a correlated multi-site treatment is neces- 
sary [53j , We do not study that scenario in this paper. 

Following Rcf. [1J] we realize that, in addition to the 
occupation number basis, cigenfunctions of the operators 
n, J-~ — J- ■ J-, and J- z also form a complete basis for the 
on-site Hilbert space of H al . In fact, these angular- 
momentum basis states \n, F, M) diagonalize 7J final when 
5 = with energy spectrum 

E & ™\§ = 0) = ^n(n - 1) + y ( F ( F + !) - 2 ") > ( 5 ) 

where n is the local atom number and the quantum num- 
ber F is restricted to F < n and even/odd F for even/odd 
n. The integer M is the magnetization quantum number 
with \M\ < F. For 8 ^ the quadratic Zeeman interac- 
tion couples the angular-momentum basis states. Find- 
ing its matrix elements is involved, leading us to perform 
all simulations in the occupation number basis. 



C. Observables 

To analyze the non-equilibrium dynamics of our sys- 
tem, we follow several observables. The first is the atom 
number per lattice site in each spin component, (fii m }, 
which can be detected either in situ [54[ or, after release 
of the atoms from the lattice, by the Stern- Gerlach sep- 
aration method where the spin states are first spatially 
separated and then detected [16[ . The second observable 
is the visibility, which is a measure of coherence of the 
wavefunction and equals the number of atoms with zero 



momentum in the spin-dependent momentum distribu- 
tion. This is the standard quantity measured after re- 
leasing atoms from the lattice and a time of flight expan- 
sion |36ll55|. Within our simulation it is given by |(ai m )| 2 
for any site i. Finally, we study the square of the in-situ 
transverse magnetization |(.Fia:)| • Transverse magneti- 
zation can be measured by Faraday rotation spectroscopy 
which allows for continuous observation of spin popula- 
tion in a BEC dm. 



III. DYNAMICS OF AN 
ANTIFERROMAGNETIC SPINOR 

A. Evolution without magnetic field 

In this section we analyze the dynamics of a longitudi- 
nal polar superfluid ground state after the optical lattice 
strength is rapidly raised. The condensate evolves un- 
der iJ final with S = for hold time t, after which one or 
more of the observables is measured. Figure QJa) shows 
typical evolution of the in-situ population of spin compo- 
nents m = ±1, as a function of hold time. Here we use 
typical numbers for 23 Na atoms - in the initial lattice 
with Uo/(zJ) = 2, U2 = +Q.04Uq, and a mean atom- 
number per site of (h) = 1.31, and in the final lattice 
Uq = 0.1hujf,U2 = +0.04£/o, where ujf is the harmonic 
frequency near the minima of the lattice potential. (An 
inifinitesimal 8 > is applied to ensure formation of the 
longitudinal polar state.) 

Initially, all atoms are in the m = state and as time 
evolves atoms begin to appear in states m — ±1 be- 
cause of spin changing collisions, and a pattern of peri- 
odic modulation emerges with a period of /1/E/2, while 
conserving zero magnetization. This spin-mixing time 
trace also contains information about the composition 
of the initial many-body state. To explore this, we fur- 
ther analyze the dynamics in Fig. UJb) and (c). Figure 
[IJb) shows the Fourier analysis of the time trace in panel 
(a), and the inset shows on-site Fock state probabilities 
|cfj| 2 for the initial state. The shape of the initial state 
number distribution is characteristic of a nearly coherent 
or a slightly-squeezed state. In the frequency spectrum, 
peaks are observed at frequencies that are integer mul- 
tiples of l/2/h. In fact, they occur at 3U2/I1, hUijh, 
and a small contribution at ITJ^Jh. These features can 
be understood from an analysis of the eigenenergies in 
Eq. flS} of H final at 6 = 0. Similar to number Fock state 
composition of the initial state, it is also a superposi- 
tion of angular momentum states \n, F, M). The ob- 
servables h m commute with total atom number h and 
thus only measure the coherence between states with dif- 
ferent F but the same n. For states with n = 2 the 
two allowed F are and 2 with energy difference 3U2- 
This leads to a peak at iU^/h in Fig. QJb). For n = 3 
states, F = 1 and F = 3 exist leading to the frequency 
at 5U2/h. The small feature at lU^jh, indicates the pres- 
ence and mixing of \n, F, M) = |4, 4, 0) and |4, 2, 0) states. 
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FIG. 1: (color online) Spin-mixing dynamics after a sudden 
increase of lattice depth starting with the longitudinal anti- 
ferromagnetic ground state of spin-1 23 Na in a zero magnetic 
field, (a) The on-site atom number of the spin components, 
(n m= o,±i), as a function of hold time showing spin- mixing. 
Time is in units of /1/C2 and system parameters are described 
in the text, (b) Fourier spectrum of the time trace in (a) show- 
ing the frequencies involved in the spin-mixing dynamics. The 
inset shows the initial ground state Fock state number proba- 
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as a function of atom number no in spin 



component m = 0. There are no atoms with m — ±1. (c) 
Contribution Po to (n m= o) of having no atoms in the m = 
state, where Po(no) = X/n _ n o|c«| 2 - The curves are for 
no = 1,2,3, and 4. It shows that the Sh/U2 periodicity is 
due to no = 2 and 5I1/U2 periodicity due to no = 1 and 3, 
respectively. 



The above analysis of the eigen energies is confirmed in 
Fig. [TJc). It depicts the time evolution of the contri- 
bution Po to (n m =o) of having no atoms in the m = 
state, where Po(no) = ^„ no|cfl| 2 . Atom number 

rig = 2 has a period of (h/U2)/3 as we oscillate between 
states \n-i,no,ni) = |0,2,0) and |1,0, 1) with a total 
of two atoms, while that for uq = 1 and 3 has a period 
of (/i/C/ 2 )/5. Here, we oscillate between the three atom 
states |n_i,n ,ni) = |0,3, 0) and |1,1,1). 

The spin-mixing dynamics can be compared and con- 
trasted with other spin-1 experiments. In Ref. [24[, a 
pair of F = 1 87 Rb atoms were prepared in a single site 
of a deep optical lattice in the Fock state |0, 2, 0), and al- 



lowed to spin- mix with |1,0, 1). Spin mixing oscillations 
between two levels analogous to Rabi oscillations were 
observed with a single frequency. On the other extreme, 
spin-mixing for a spinor BEC with a large number of 
atoms has been discussed in theory and observed in ex- 
periments [lj, LL5| ■ They are in a regime where a classical 
pendulum phase-space analysis is appropriate [571 . |58| . 
and although there can be spin-mixing oscillations for 
specific initial states, quantum recurrences due to the 
discrete energy spectrum is absent. Our analysis here ex- 
plores the regime which is between these two - the single 
Fock state regime and the regime of large atom number 
condensate. As such, we are exploring a regime which 
can shed light on the semi-classical transition to large 
condensate dynamics, a topic for future investigation. In 
our analysis here, we can analyze the multiple frequencies 
of the dynamics time trace to probe the composition and 
atom number statistics of the initial many-body state. 

Figure [2] shows the dynamics of visibility of the m = 
state - the occupation of the zero momentum state for 
m = component, for the same initial state and param- 
eters as in Fig. [TJ The visibility \(a m= o(t))\ 2 measures 
the phase coherence in this spinor superfluid system. 
We show the relative visibility |(a m= o(£))| 2 /|(a m= n(£ = 
0)}| 2 in Fig. [2ta). We see that the atoms oscillate be- 
tween being completely coherent to completely incoher- 
ent (|(a m =o)| 2 ~ 0). The pattern is more complex than 
the spinless boson visibility |35|. There is a fast oscil- 
lation with time scale ~ h/Uo, which is modified by a 
slower envelope with a time scale ~ hjXJ?,. The exact 
nature of the complex oscillations is revealed in the fre- 
quency spectrum shown in Fig. &h). Similar to Fig. [TJ 
features appear at small integer multiples of TJijh. Here 
they are located at 2U2/I1, "SU^/h, and 5Uz/h. Peaks 
also occur at much larger frequencies with a dominant 
frequency at (Uq + U2)/h, which is 26U2/h in this ex- 
ample. Twenty six is also the number of fast oscillations 
in a full period hjU-x as can be seen in panel (a). This 
indicates that for an unknown ratio U2/U0, one full re- 
vival oscillations of the visibility can help determine this 
ratio by counting the number of fast oscillations or equiv- 
alcntly performing a frequency analysis. Combined with 
the realization that this ratio is independent of lattice pa- 
rameters, this method of determining spinor interactions 
is one of the key findings of this paper. 

The visibility spectrum frequencies appear because the 
expectation value of the annihilation operator is sensitive 
to the overlap between Fock states of different atom num- 
ber [3g. For example |n_i,nn,ni) = |0,2,0) connects to 
1 0,1,0) giving rise to the dominant frequency £/n + U2. 
Other frequencies can similarly be explained by perform- 
ing an expansion of the initial state in the angular mo- 
mentum basis, and applying the evolution operator for 
the final Hamiltonian Eq. (TTJ) at S = 0. For higher occu- 
pation numbers not shown here we find that the visibility 
patterns become more complex. In fact, by controlling 
the initial squeezing [36l . |52| (i.e. by controlling the ini- 
tial tunneling energy J) as well as the total occupation, 
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FIG. 2: (color online) (a) Visibility |(a m= o(t))| 2 /|(a m=0 (i = 
0))| 2 of the m — state as a function of hold time t after a 
sudden increase of the lattice depth, for the same parameters 
as in Fig. [TJ The pattern of collapse and revival of coherence is 
more complex than the spinless boson case due to the compe- 
tition between spin-dependent and spin-independent interac- 
tions. There are U0/U2 + I fast oscillations in one full collapse 
and revival period of h/Ui, yielding a method to determine 
the ratio U2/U0 from visibility oscillations, (b) Spectrum of 
the visibility oscillations showing the range of contributing 
frequencies, the most dominant one being at U0/U2 + 1 = 26. 



we can change the complexity of the frequency spectrum 
and thereby make it amenable for analysis. 



B. Effects of magnetic field 

The energies of three component spin-1 atoms are 
sensitive to external magnetic fields. Such external 
fields have been exploited in manipulating spinor atoms 
in an optical trap - to access ground state proper- 
ties [H 0| i m detection such as in Stern Gerlach separa- 
tion method 11611 . or to influence the dynamics in quench 
experiments [39(. For spinor atoms in an optical lat- 
tice external magnetic fields cannot be ignored and, in 
fact, lead to unique physics. For example, the quadratic 
Zeeman shift affects the phase diagram [33|]. The ra- 
tio of quadratic shift to the spin-dependent interaction 
strength, 6/U2, controls the physics of this system. Two 
different regimes emerge - the Zeeman reg ime for 8 > U2 
and the interaction regime for S < U2 |12| . 

Figure [3] shows an analysis of spin-mixing oscillations 
(n m —o) in the presence of a quadratic Zeeman shift 5, 
during the initial state preparation and during the evolu- 
tion. The other parameters are as in Fig. [TJ To highlight 
the effects of a magnetic field, we show a comparison of 
the dynamics for 5/U~2 = 0, 0.5, 1.5, 2.5, and 5 in Fig.[3J]a). 
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FIG. 3: (color online) Spin-mixing oscillations in the presence 
of a magnetic field in the form of a quadratic Zeeman shift. 
Panel (a) shows the population in m = as a function of 
hold time for several values of the quadratic Zeeman shift. 
Panel (b) shows the frequency analysis of panel (a) with two 
dominant frequencies in the time evolution. At S = these 
frequencies are 3C/2//1 and 5U2/h. The inset shows that the 
two dominant frequencies as a function of S first dips before 
slowly increasing due to the competing nature of 5 and XJi- 
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FIG. 4: (color online) (a) Spin-mixing dynamics (n m ) for a 
ferromagnetic Uz < rubidium condensate as a function of 
hold time, (b) Fourier spectrum showing the relevant frequen- 
cies of 3U2/h, 51/2/h and 7U2/h, and their relative ampli- 
tudes, (c) Visibility |(a m= n(i))| 2 /|(a m= o(£ = 0))| 2 as a func- 
tion of hold time, (d) Frequency spectrum of the visibility 
showing features at (Uo — U2)/h and 2([/n — U2)/h, pointing 
a way to measure interaction ratios U2/U0 from the visibility 
spectrum. 



three-atom Fock states. In a condensate with large par- 
ticle numbers the competition between interaction and 
Zeeman energy gives rise to a sharp phase boundary at 
6=1/2 which can be manifested through magnetic field 
induced resonances I59J. 



IV. DYNAMICS OF A FERROMAGNETIC 
SPINOR 

A. Evolution without magnetic field 



We see that the oscillations become faster while simul- 
taneously the amplitudes get smaller for increasing 5 in 
the Zeeman regime 6 > Ui- The spin-mixing dynamics 
vanishes for a large enough i?-field. In panel (b) we plot 
their frequency spectra comparing the frequencies and 
the amplitudes. The inset shows two dominant frequen- 
cies as a function of the Zeeman strength. A closer look 
at panel (b) reveals that for a S in the interaction regime 
6 < U2, the dominant frequency initially decreases before 
starting to increase, due to a competition between the 
spin-dependent interaction and the Zeeman term. Much 
of the S-field effects can be understood from the eigen- 



24] 



values of H Rnal calculated in the Fock state basis 
using conservation of atom number and magnetization. 
For up to three atoms per lattice site this involves diag- 
onalizing at most a 2 x 2 matrix. The energy splitting 
of the eigenvalues of the 2x2 matrices is observed. In 
fact, from top to bottom the two curves in the inset of 
panel (b) are due to contributions from two-atom and 



Figure 0] shows the dynamics when our initial state is 
a ferromagnetic supcrfluid state of 87 Rb with U2 < 
created in a shallow lattice. For a ferromagnetic state, 
the collective spin configuration is such that the spin- 
dependent interaction energy is maximized. This means 
that the ground state is a superposition of magnetiza- 
tion states in the angular momentum basis \n,F,M) 
with variance Am ^ although it still has magne- 
tization (m) = 0. The parameters are (n) = 1.84, 
U /(zJ) = 2 and U 2 = -0.005L/ in the shallow lat- 
tice and Uo = O.lHuif and U2 = — 0.005[/n in the final 
deep lattice. The population dynamics (n m ) is shown 
in Fig. QJa). The oscillation amplitude is not as large 
as in the polar case in the previous section. There are 
two reasons for this: first, the ground state has compa- 
rable populations in all three components and therefore, 
the population difference between the spin components 
is smaller to begin with, unlike for 23 Na where initially 
only the m = state is occupied. Second, the initial 
state is much closer to the ground state in the deep lat- 



tice. We would like to point out that the smallness of 
the ferromagnetic interaction is not responsible for the 
oscillation amplitudes being small. A Fourier analysis 
in panel (b) of the m = population shows that the 
frequencies present are 3U2/h, §U-}/h and lUz/h, as in 
the polar state. Again, energy differences obtained from 
Eq. ([5]) give us those frequencies. These frequencies and 
their spectral weight determine the composition of differ- 
ent Fock components in the initial many-body state, and 
therefore, can be used as an experimental probe. We will 
show in the next subsection that adding magnetic-fields 
can enhance the amplitude of spin-mixing dynamics. 

The visibility \{a m =o(t)}\ 2 /\{a m =o(t — 0))| 2 is shown 
in Fig. |4|c). It has a simpler pattern than for the po- 
lar case in Fig. [5] The coherence of the initial matter- 
wave exhibits collapse and revival modulations with a 
fast timescale of h/Uo- Its frequency spectrum in panel 
(d) shows two dominant frequencies, {Uq — V^jh and 
2(£/n — U2)/h, although other frequencies with extremely 
small amplitudes do exist. For 87 Rb the dominant fre- 
quencies are 199 and 398 in units oiU^jh. As explained 
for 23 Na, these frequencies appear from the overlap of 
Fock states connected by the annihilation operator a [36[ . 
Here, the peak at (Uo — U2)/h appears due to the overlap 
of number states |0, 1, 0) and |0, 2, 0) . Similarly, the peak 
at 2(Uq — U2)/h appears due to number states |0, 2, 0) and 
|0, 3, 0). As with 23 Na, finding these frequencies yields a 
method to experimentally determine VijU^. Neverthe- 
less, since we need to observe many oscillations, « 200 
for 87 Rb and ~ 80 for 41 K, this will be a challenging 
application of quantum phase revival spectroscopy. 



B. Effects of magnetic field 

In this subsection we study the effects of the quadratic 
Zeeman interaction on a ferromagnetic spinor. In Fig. [5] 
we show the dynamics when the quadratic Zeeman shift 
is S = U2 and I.5E/2 for the top and bottom row re- 
spectively. We plot spin population (n m =n) an d in-situ 
transverse magnetization |(.Fa;)| 2 for 87 Rb. The parame- 
ters are (n) = 1.31, U /(zJ) = 2 and U 2 = -0.005*7 in 
the shallow lattice and Uq = O.lfiw/ and U2 = — 0.005E/n 
in the final deep lattice. In all the cases, the oscillations 
are no longer periodic in hjU^-, but the values of 5, Uq and 
U2 combined with the initial state composition influence 
the dynamics. We find that the spin-mixing modulation 
amplitudes become large in the presence of a B-ficld un- 
like the antiferromagnetic case in the previous section. 
The oscillations also get faster, and the fast spin oscilla- 
tions are modified by an envelope of a slower modulation 
pattern which is manifested in a beat-like pattern involv- 
ing the two dominant frequencies as clearly evident in 
Fig.[5jc). As we increase S, the population differences in 
the m = and m = ±1 states become larger, and there- 
fore a large spin-mixing modulations can occur. Trans- 
verse magnetization |(.F X )| 2 shown in (b) and (d) goes 
through partial revivals and complete collapse. Faraday 
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FIG. 5: (color online) Effect of quadratic Zeeman interactions 
for 87 Rb dynamics for 8 = 1/2 and I.5U2 for the top and bot- 
tom row respectively, (a) and (c) shows spin mixing dynam- 
ics of the m = state, (b) and (d) shows in-situ transverse 
magnetization \(T X }\ 2 oscillations. Spin-mixing and K^)! 2 
oscillations depend on a combination of 8, Uq, U2 and the 
initial state. We see that the spin-mixing amplitudes increase 
as we increase the field. On the other hand transverse mag- 
netization oscillations get weaker. 



rotation spectroscopy can be used to detect transverse 
magnetizations Q. This could give us a direct probe of 
the magnetic properties in our quench set up, in addition 
to the coherence properties which we can obtain through 
population and visibility revivals. We show in panels (b) 
and (d) that the initial amplitude for \{F X )\ 2 decreases for 
increasing B-field, but the number of collapse sequences 
increases. When 5 > 2C/ 2 the transverse ferromagnetic 
supcrfluid turns into a longitudinal superfluid for which 
|(J r a; )| 2 =0. In that regime transverse magnetization is 
zero throughout the evolution, but large amplitude spin- 
mixing modulations still occur. 



V. SIGNATURE OF EFFECTIVE THREE-BODY 
INTERACTIONS 

A. Three-body interactions 

In our dynamical scheme, the evolution of the many- 
body state takes place in a deep optical lattice with 
no tunnelling to the neighbors. In such a setting, 
even within the single-band Bosc-Hubbard model, there 
are effective three and higher-body interactions due to 
collision-induced virtual excitations to higher bands or 
vibrational levels. For the spinless bosonic case, such 
effective multi-body interactions have been predicted in 
theory [371 ]. A recent experiment [36( monitored matter 
wave collapse and revival dynamics for tens of oscillations 
and observed the signature of higher-body effects in the 
visibility time trace. A more accurate treatment of a 
multi-component system in a deep lattice should there- 
fore also incorporate higher-band induced multi-body in- 



teraction terms. Three-body interactions can be impor- 
tant for Efimov physics |60j and for many-particle sys- 
tems giving rise to novel and exotic phenomena [6lT - l63j . 
In a deep lattice, the minimum of the potential at a 
single site can be approximated as a harmonic potential. 
For spin-1 bosons in a single isotropic harmonic trap, the 
derivation of the effective three-body interactions is given 
in the appendix. The effect can be concisely represented 
by adding 



H 3B .ctt = -7-n (n 




l)(n-2) + ^-(f 2 -2n)(h-2) (6) 



to H al in Eq. |4]). Here, as previously, n = fi-i+no+fii 

is the on-site atom number and _F = J2 m m' a \rJ ? rnm l <W 
is the on-site total angular momentum. The Hamilto- 
nian term with strength Vq is similar to the spin-0 ef- 
fective three-body interaction of Refs. |36|, [37| and only 
depends on total particle number. The new Hamiltonian 
term with strength V 2 depends intricately on the atomic 
spin. For a harmonic potential with frequency w*, Vq is 
attractive and equal to V = -1.34C/^ /(hujf) [37]]. The 
strength V 2 satisfies V 2 = 2(U 2 /Uo)V~o- There is no mag- 
netic field dependence in the effective three-body terms. 
Finally, we note that the perturbativc effective three- 
body Hamiltonian is only valid when the three-body in- 
teraction strengths (Vo,V 2 ) are much smaller than the 
corresponding two-body strengths (Uq, U 2 ). 

The on-site effective Hamiltonian i?3B, c ft is diagonal in 
the angular momentum basis states \n,F,M) with diag- 
onal matrix elements 



E 



3B,cff 



Vb 

G 



n(n-l)(n-2) + — (F(F + l)-2n)(n- 



•2). 

(7) 

with the same restrictions as in Eq. ([5]) on allowed values 
of F. Equation (0 is one of the key results in this paper. 
This spectrum extends and generalizes the spectrum of 
spin-mixing spinor hamiltonian as presented in Ref. [I4j | 
to the effective three-body case. As discussed next, this 
helps us quantify and understand the effects of three- 
body interactions on our dynamics scenario. 



B. Effects on Dynamics 

For spinless bosons in an optical lattice, a time of flight 
measurement of the visibility dynamics determined the 
strength of the effective multi-body interactions [3|| |52| ■ 
Here for spin-I bosons in an optical lattice, we show 
that the effective three-body interaction effects can be 
observed directly in the on-site population density - in 
the spin oscillation dynamics. This opens up the possi- 
bility that in situ measurements such as quantum micro- 
scopes [5J, |64j and other techniques |65[ in lattices could 
be used to detect effective multi-body interactions. Be- 
cause of the more complex nature of visibility patterns, 
we only analyze spin-mixing population dynamics. 
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FIG. 6: Effects of effective three-body interactions on the 
spin-mixing dynamics for a polar ( 23 Na) initial state, (a) The 
population (ft m= o) as a function of hold time for occupation 
(n) = 2.35 and 5 = with (solid line) and without (dashed 
line) the effective three-body interaction, (b) Frequency anal- 
ysis of the time trace with effective three body interactions 
reveals the presence of additional frequencies, which could be 
used to determine the three-body interaction strength. 



In Fig. [6] we show the effects of the effective three-body 
interactions in the spin-mixing dynamics of an initial po- 
lar superfluid state of 23 Na at U /(zJ) = 2, U 2 = 0.04C/ 
and 6 = 0. The occupation is (h) = 2.35, which high- 
lights the three-body effects, and (m) = 0. The interac- 
tion strengths for the deep lattice are Uq = O.lhujf and 
U 2 = 0.04(7 , so that V = -0.134C/ and V 2 = -0.268C/ 2 - 
In an earlier section we have seen that spin-mixing dy- 
namics is controlled by U 2 - Here, we find that this spin- 
mixing scaling is also influenced by the three-body in- 
teraction V 2 . Figure Ufa) shows a comparison of the dy- 
namics with and without the three-body term. We see 
that the time traces start to differ after the first oscil- 
lation, and the periodicity in h/U 2 is destroyed. A fre- 
quency analysis of the oscillations in panel (b) elucidates 
the exact nature of the modulations. Without a three- 
body term, strong features appear at 3U 2 /h, 5U 2 /h, and 
7U 2 /h. With a three-body term, additional frequencies 
appear at 2A6U 2 /h,4.55U 2 /h, and 5.75U 2 /h, which fol- 
low from the energy differences of the three-body spec- 
trum in Eq. f7j). Identification of any of the frequencies 
gives us the V 2 coupling strength. For example, the peak 
at A.55U 2 /h arises as the initial state contains contribu- 
tions of angular momentum states |3, 1,0) and |3, 3,0) 
containing three atoms. These two states have an energy 
difference of 5U 2 + 5V 2 /3 whose signature is the 4.55U 2 /h 
frequency. 

In experiments an unknown effective three-body 
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FIG. 7: Effects of effective three-body interactions on the 
spin-mixing dynamics for a ferromagnetic ( 87 Rb) initial state, 
(a) The population (n m= o) as a function of hold time with 
(solid line) and without (dashed line) the three-body inter- 
action. For 87 Rb the oscillation amplitude is small, (b) 
Frequency analysis reveals the presence of several frequen- 
cies, which are a signature of induced three-body interactions; 
compared to 23 Na in Fig. [5] one frequency is missing for 87 Rb 
due to its ferromagnetic nature. 



strength V*2 can be deduced by assigning several frequen- 
cies in the time trace. The presence of all the other fre- 
quencies can be used to reduce error in the measurement 
and verify the spectrum Eq. ((?]). The spectral weights 
and the values of the frequencies also give us clue about 
the initial superfluid state. Control of initial squeez- 
ing [3a, [52|, by varying Uq/{zJ) in the shallow lattice, 
can be used in such a way that some of the frequencies 
are more dominant so to make detection easier. 

In Fig. [7] we show the spin dynamics (h m —o) for a fer- 
romagnetic ( 87 Rb) initial state at Uo/(zJ) = 2, 6 = 0, 
(rh) = 0, and (h) = 2.35 in the shallow lattice. In the 
deep lattice we use Uq = O.lhajf and U2 = — 0.005l7 . 
The spin-mixing amplitude is not that prominent for 
87 Rb. Nevertheless, the influence of the three-body 
interaction is discernible in our simulation. The fre- 
quency spectrum of the time trace shown in (b) makes 
it clear that new frequencies emerge at 4.55E/2//1 and 
5.75U2/h proving the existence of three-body interactions 
and yielding a method to measure its strength. The ap- 
pearance of additional frequencies is similar to the Na 
case, except that we do not observe a feature at 2A6U2/I1, 
which is due to a coherence between four-atom angu- 
lar momentum states |4, 0, 0) and |4, 2, 0). For the ferro- 
magnetic initial superfluid the angular momentum state 
|4, 0, 0) is absent. It is conceivable that initial squeezing 
control of a ground state [3|| [52| or other specific initial 
state preparations [lj, [l5( can be used to see larger am- 
plitude spin-mixing oscillations to make the experimental 



detection of three-body effects easier for ferromagnetic 
coupling. 



VI. CONCLUSION 

In this paper we performed a theoretical study of 
the dynamics of spin-1 bosons in an optical lattice in 
a quench scenario where we start from a ground state 
in a shallow lattice and suddenly raise the lattice depth. 
We have shown that the ensuing spin-mixing and visi- 
bility oscillations can be used as a probe of the initial 
superfluid ground state. The spectral analysis of time 
evolution reveals the Fock state composition of the ini- 
tial state and thereby its superfluid and magnetic prop- 
erties. Analysis of visibility oscillations, i.e. quantum 
phase revival spectroscopy, further yields a method to 
determine the spin-dependent and spin-independent in- 
teraction ratio U2/U0, which is an important quantity 
for spinor gases. We treat both antiferromagnetic (e.g. 
23 Na atoms) and ferromagnetic (e.g. 87 Rb) condensates. 
For ferromagnetic interactions the spin-mixing oscillation 
amplitudes are small. When external magnetic field can- 
not be ignored, the inclusion of a quadratic Zeeman field 
is necessary, and we have quantified such dynamics. We 
have shown that the presence of a magnetic field increases 
the spin-mixing amplitude for a ferromagnetic conden- 
sate. 

The Hamiltonian that more accurately describes the 
physics of the final deep lattice comprises of two-body 
as well as effective multi-body interactions, which arise 
due to virtual excitations to higher bands. We derive 
the induced three-body interaction parameters for spin-1 
atoms in a deep harmonic well and show its effect on the 
spin-mixing dynamics. We demonstrate that a frequency 
analysis of the oscillations can detect the signature and 
strength of the spin-dependent three-body interactions. 
We stress our finding that the three-body interactions 
for spinor atoms can be observed directly in the in-situ 
number densities in addition to the time of flight visibility 
as observed for spinless bosons |36j . 

Although there have been many theoretical studies for 
spin-1 bosons in an optical lattice, a many-body corre- 
lated ground state has not yet been achieved experimen- 
tally. There are many unexplored questions in that re- 
gard. Here, we have combined the study of dynamics 
with optical lattice spinors to show how non-equilibrium 
dynamics can be used as a probe for revealing ground 
state properties and spinor interactions. There are other 
dynamic scenarios that can give different perspectives on 
spinor lattice physics, such as a quench from Mott insu- 
lator to superfluid, evolution in a tunnel coupled lattice, 
to name a few. Interplay of superfluidity, magnetism and 
strong correlations makes this a rich system where study 
of quantum dynamics may lead to a better understanding 
of collective phenomena. 



10 



Appendix A: Effective three-body spinor 
interactions 



In this appendix wc derive the effective three-body 
spinor interaction due to virtuaf excitations to excited 
trap ieveis of the isolated sites of the dee p op tical lattice. 
The derivation closely follows Refs. [33, [52(. Atoms are 
held in the ground state of a site, which we approximate 
by an isotropic 3D harmonic-oscillator potential with fre- 
quency ujf. Then the Hamiltonian for spin-1 bosons is 
H = H0 + V0 + V2 with 



H = (e a + m 2 S)al na a ma (Al) 

Vb = 7; c ogaPr<s a ka a w a kaaii3 (A2) 

V2 = 7;C2gaPr<S a ka a w( F krn ■ Fi n )a m ~fa n s , (A3) 



where operators a ma annihilate an atom in spin projec- 
tion m = —1,0,1 and 3D harmonic oscillator state a. 
Wc use the convention that roman and greek subscripts 
represent spin projections and harmonic oscillator states, 
respectively. Repeated indices are summed over. The 
single-particle energies e a and m 2 S are the harmonic os- 
cillator and quadratic Zeeman energies, respectively. For 
the ground state a = 0. We choose eo = and will re- 
quire that e a — eo -C S for a =/= 0. The spin-independent 
and spin-dependent atom-atom interactions Vo and V2 
have "bare-coupling" strengths cq and C2, respectively. 
The vector F = (F x ,F y ,F z ) are spin-1 matrices. Each 
symmetric real coefficient g a p rf s is a 3D integral over the 
product of the four oscillator wavefunctions a, j3, 7, and 
5. In this appendix the two atom-atom interactions are 
explicitly normal-ordered in order to facilitate the deriva- 
tion. 

We can now derive the effective Hamiltonian with 
atoms in the lowest oscillator state using degenerate per- 
turbation theory with zcroth-order Hamiltonian Hq and 
perturbation Vq + V2. First, wc define for every spa- 
tial mode a the spin wavefunction \{n—i,riQ,ni} a ) with 
n m = 0,1,2,... atoms in spin state m. The ground 
states P are formed by the orthonormal basis functions 
\g) = |{ji_i,no,rii}o,0 a ^o) f° r anv value of n m and 
where a ^o indicates that there are no atoms in excited 
spatial modes. Their energy is E g = (n_i + n\)8. Ex- 
cited states |e) with energy E e are states where at least 
one atom occupies aa/0 spatial mode. 

We reproduce the Hamiltonian in Eq. ^ in first-order 
perturbation theory once we make the assignment Uq = 
coffoo,oo and U2 = C25oo.oo- To second-order in degenerate 
perturbation theory the matrix element for ground states 
Iff') and \g) is 



(g'\SH^\g) = lj2(9'\ V o + V 2 \ 



e^P 



1 



Eg - E e 



(A4) 



Eg' - E e 



(e\V + V 2 \g) 



The sum over excited states can be evaluated by realiz- 
ing that only states |e) ex a^ a aLa m oa„o|s') with a/3 7^ 00 
and k + l = m + n contribute as both Vo and V2 conserve 
atom number and magnetization, and can only change 
the state of two atoms at the same time, i.e. states |e) 
are those with at most two atoms in the higher trap levels. 
By inspection, we then realize that the energy differences 
E e — Eg sa e a + ep are independent of the total number 
of atoms in the ground states and to good approximation 
are also independent of the quadratic Zeeman energy, as 
5 <C e Q for a^0. Similar expressions hold for E e — E g i. 

Inserting the (normalized) expression for |e) and per- 
forming the sums over |e) as well as those appearing in 
the potentials Vq and V2 we first find 



W\5H^\g) 
1 






1 



ep 



-gap, 00 



x \ j%{g^ a rnO a nO a ^ a nual a a}pa m aio\g) 



(A5) 



1 



CoC2{g'\a} p0 a\ {F pm ■ Fo^am^an^a^ajpakoaiolg) 



where the remaining sums over trap levels have been 
made explicit, repeated roman indices are still summed 
over, and we have omitted the contribution propor- 
tional to V2 times V2 as the spin-dependent interaction 
strength is an order of magnitude smaller than the spin- 
independent one. 

By normal ordering the creation and annihilation oper- 
ators in Eq. (|A5|) and using that ground states P contain 
no atoms in excited trap levels we find 



g'\6HW\g) = (g'\SH 2B + H 3B \g) 



(A6) 



where 
SH 2B 



Z. 



■{■ 



; c o a lo a lo a koaio 



(A7) 



-c C2al n0 a'l (F mk ■ F n i)a kQ ai > 



is a correction to the pair-wise two-body interaction and 
\ c o a Io a K) a mo a fcoa2oa m o (A8) 



-2c C2al al n0 a l n0 (F mk ■ F n i)a ko ai a o0 > 



is an effective three-body interaction. Here 

1 



Zi 



E 

M^OO 



5oo „ 



+ e u 



gnu,oo 



and 



2^ff00,/i0 — SpO.OO ■ 

^0 ^ 



(A9) 



(A10) 



The sums in Zn diverge and must be regularized and 
renormalized [52| . |66| . That is we require that the bare 
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coupling constant cq is denned such that, by combining 
the first- and second-order contributions, Co3oo,oo — c o-^2 
is finite and equal to Uq- Similarly, we require that 
C2ffoo,oo — 2cqc 2 Z 2 is finite and equal to U 2 . 

The sums in coefficient Z3 of the effective three-body 
Hamiltonian do converge. Hence, we redefine the effec- 
tive three-body interaction as 



1, 



•t „t J 



Hsb — -Vq a ko a lo a mQ akoaioa m0 



(All) 



+ 7T^2 a oO a mO a no(Fmk ' F n i)a,koaioa , 

where, consistent within our second-order pcrturbativc 
calculation, 

V = -6C/ 2 ~^—Z 3 and V 2 = -12U U 2 Z 3 = 2^V 
3oo,oo ^0 

are the spin-independent and spin-dependent three-body 
interaction strength, respectively. For a spherically sym- 
metric harmonic oscillator Vq = —1.34 . . . £/q /(&*;/) [371 ]. 
Since U 2 > for Na atoms, V 2 < 0, while for 87 Rb U 2 < 
and thus V 2 > 0. 

The two- and three-body interaction can be rewrit- 
ten in terms of the angular momentum operators J- = 



(P X ,T V ,T Z ) defined in Eq. flTJ) following Ref. [ij]. By 
combining the quadratic Zeeman interaction as well as 
the two-body and effective three-body interaction we find 
our final result 



H cfi = 5^m 2 a\ n a m + -U h{h-l) + -U 2 {F 2 -2h) 

m 

+ ^V h(n - l)(n - 2) + \v 2 {P - 2n)(n - 2) , 



where we suppressed the ground-state oscillator index, 
and T 2 and the number operator h commute. 
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